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Abstract. A method to evaluate spin networks for (2+l)-dimensional quan- 
tum gravity is given. We analyse the evaluation of spin networks for Lorentzian, 
Euclidean and a new limiting case of Newtonian quantum gravity. Particular 
attention is paid to the tetrahedron and to the study of its asymptotics. More- 
over, we propose that all this technique can be extended to spin networks for 
quantum gravity in any dimension. 

1 Introduction 

In this paper we study the evaluation of spin networks in (2+1) -dimensional 
Lorentzian, Euclidean and Newtonian quantum gravity. The concept of spin 
networks was first introduced by Penrose in jTT], in order to describe a combi- 
natorial picture of the geometry of space-time. 

The evaluations of spin networks are of great importance in the framework of 
Spin Foam models of quantum gravity. The first example of a Spin Foam model 
was proposed exactly in (2+1) dimensions. This is known as the Ponzano-Regge 
model ^2], and it is a model of Euclidean quantum gravity. This is constructed 
in a triangulated (2+l)-dimensional manifold M, by assigning an irreducible 
representation of 577(2) to each edge of the triangulation. The corresponding 
representations are labelled by half integers j = 0, |, 1, .... Then to each tetra- 
hedron we assign a 6j-symbol. We then multiply all of the 6j-symbols assigned 
to the tetrahedra of the triangulation of our manifold and then sums over all 
the irreduclible representations of our group SU{2). 
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This is a kind of discretised path integral for Euclidean quantum gravity. 
Ponzano and Regge then studied the asymptotics of this Euclidean 6j-symbol. 

One question one can ask first is whether a similar Ponzano-Regge formula- 
tion exists for (2+l)-dimensional Lorentzian quantum gravity and whether we 
can study the asymptotics of the Lorentzian tetrahedron network. This problem 
was studied first in jH] and in [7], by analysing the asymptotics of a Lorentzian 
tetrahedron labelled by discrete unitary representations of the Lorentz group. 
Then in jH] a state sum model for (2+l)-dimensional Lorentzian quantum grav- 
ity was proposed. Later the asymptotics of the square of the Lorentzian 6j- 
symbol was studied in P] . One of the new problems we study in this paper is the 
asymptotics of the usual Lorentzian tetrahedron network labelled by continuous 
unitary irreducible representations of the connected Lorentz group 5*00(2,1). 
Inspired by |3j we study, in a similar fashion, Lorentzian (2+l)-dimensional 
quantum gravity. Particularly we pay attention to the tetrahedron which is of 
great importance in the evaluation of transition amplitudes and of the parti- 
tion function. There are different kinds of tetrahedra in Lorentzian geometry 
according to whether their corresponding edge vectors are timelike, spacelike or 
null. When dealing with spin foam models for quantum gravity the language of 
representation theory comes into play and the Lorentzian tetrahedra spin net- 
works may have edge vectors labelled by any mixture of discrete, or continuous 
representations of the connected Lorentz group SOq(2,1). When considering 
the continuous principal unitary representations only we can attempt to con- 
struct a spin foam model for (2+l)-dimensional Lorentzian quantum gravity by 
using the dual complex J7a to a triangulation A of our 3-dimensional manifold. 

In order to construct the model we label each face of the complex j7a by 
a principal unitary irreducible representation of the connected Lorentz group 
SO (2, 1), or equivalently if we think of the triangulation, we label each edge 
instead. The state sum model is then given by 



where the integration is carried over the labels of all internal faces of the dual 
complex and A(f), A(e), A{v) are the face, edge, and vertex amplitudes respec- 
tively. 

These amplitudes are given by the evaluation of spin networks such as 




(1) 




A(e) 
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A(v) 




p 



34 



The definition for evaluating spin networks for this (2+l)-dimensional Loren- 
ztian model is one of the problems which concerns us in this paper. In particular 
it is the evaluation of the tetrahedron network (vertex amplitude), that we con- 
centrate on. 

The way to evaluate our spin networks is based on a similar idea first de- 
scribed in for the evaluation of relativistic spin networks. This method was 
then applied to the study of the 4-simplex in |5J. The idea was interpreted to 
the case of simple spin networks in ^0] by thinking of them as Feynman inte- 
grals over an internal space. This idea is given by using propagators (or kernels) 
which are defined in the internal space. Then with every edge of the spin net- 
work we associate a propagator, and then we take a product of all these kernels 
and integrate over a copy of the internal space for each vertex. 

Particularly for our Lorentzian case model the above recipe, applied here 
for the first time, translates as follows: The first ingredient we need for our 
description is a kernel K p (x,y) where p is a real number which stands for a 
unitary principal irreducible representation of SOq(2, 1) and x and y are two 
points in the hyperbolic plane H 2 . The latter can be seen as the upper sheet of 
a timelike two-sheeted hyperboloid in 3-dimensional Minkowski space-time. 

Once we have this kernel K p (x,y), in particular the evaluation of the tetra- 
hedron will be given by the following integral 



where xi, x 2 , x 3 , x 4 G H 2 and our integral is computed over 3 copies of the 
hyperbolic plane 1 . 

One of the problems is to know whether our tetrahedron evaluation con- 
verges. Once this is true our next step is to study the asymptotics of the above 
integral. 

This way of evaluating spin networks gives us an attempt to define the 
computation of transition amplitudes between spin networks by using the spin 

1 Originally we would have an integral over 4 copies of H 2 , but regularisation requires that 
one of these integrals is dropped 




dx 2 dx- i dx A K PX2 (x 1 , x 2 )K pi3 (xi, x 3 )K Pl4 (x 1 , x 4 ) 



xK p23 (x 2 , x 3 )K p24 (x 2 , x 4 )K P34 (x 3 , x 4 ) 



(2) 
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foam model description of formula (1). This gives us the analogue of the previous 
known Euclidean Ponzano-Regge model [T2"] . 

Then it is immediately logical to ask whether we can apply the same tech- 
nique to the Euclidean model and to the asymptotics of the tetrahedron network. 

This is then our next step, to study in a similar fashion the case of Euclidean 
quantum gravity whose group is given by SO (3). The irreducible representa- 
tions of SO (3) are given by positive integers. Then again, if we think of the dual 
complex of our triangulation, we can define a spin foam formula for Euclidean 
quantum gravity given by 

oo 

z(m) =y,u A (f) n a ^ n a m (3) 

j=0 f e v 

This is the analogue of our integral formula (1), and moreover, it is very similar 
to the Ponzano-Regge formula. The important difference is to notice that in the 
Ponzano-Regge model we have the group £77(2) which is the double covering 
of £0(3). Then we go on to calculating the amplitudes for Euclidean spin 
networks by using the described recipe. This case has been studied before in 
|10j . However our treatment completes the whole story, as we include all the 
details of the calculation of the asymptotics of the tetrahedron network by using 
the stationary phase method and moreover, we mention the case of degenerate 
configurations. 

The evaluation goes as follows: considering a kernel kt(x, y), where £ denotes 
an irreducible representation of £0(3). For this case x, y are points on the two 
dimensional sphere £ 2 . Our tetrahedron network will be given by the integral 



/ dx 1 dx 2 dx 3 dx 4 Ke 12 (x 1 , x 2 )K £l3 (x 1 , x 3 )K eii (xi, x 4 ) 

7(5 2 ) 4 

XKl 23 (x 2 , X 3 )Ki 2i (x 2 , £ 4 )7Q 34 (£3, x 4 ) (4) 

In this case we do have an integral over 4 copies of £ 2 and no regularisation is 
required, as £ 2 is compact. 

We complete our picture by studying a completely new case, at least in the 
(2+l)-dimensional quantum gravity case, which we call Newtonian quantum 
gravity. This case deals with IS 0(2) spin networks. The construction of such 
kind of quantum gravity is analogous to the Lorentzian principal unitary 

irreducible representations of IS 0(2) are also labelled by positive real numbers 
M + . The evaluation of the tetrahedron network for this case is given by a similar 
integral to integral (2), with the only difference that for the Newtonian case, 
we integrate over three copies of the flat plane R 2 . 

The idea is that all these three cases are related to each-other by a limit 
process, and then we have a whole picture of a unified theory. This unified 
theory idea was also proposed in pQ for the (3+l)-dimensional quantum gravity 
case. 
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The main questions we ask then are: What do the kernels look like? What 
are the evaluations of the tetrahedron networks in each case? 

We divide this paper in the following sections. In section 2 we study the case 
of (2+l)-dimensional Lorentzian quantum gravity spin networks. The compu- 
tation of the spin networks with the help of an internal space is completely new. 
The whole idea and techniques were discovered by following a similar method 
to the one used for the (3+1) -dimensional Lorentzian Barrett-Crane model. 

In section 3 we study with the same spirit the model of Euclidean spin 
networks. 

In section 4 we study a new kind of model which we call Newtonian spin 
networks. The kind of spin netwoks for this model were studied for the case of 
4-dimensions in pQ. However for (2+l)-dimensions is completely new as well as 
the name we give to it. The explanation of why we call them Newtonian is also 
new. 

Finally in section 5 we conclude by proposing a way to study spin networks 
for quantum gravity in any dimension following the spirit of the paper and which 
generalises our (2+l)-dimensional case. Furthermore for (3+l)-dimensions we 
just have again the well known Barrett-Crane evaluations. 

2 Lorentzian spin networks 

Equation (1) gives a way to define a partition function for the (2+l)-dimensional 
Lorentzian quantum gravity as a spin foam model, for the case in which we are 
only considering the principal unitary representations of the group 5*00(2, 1). 

In this section we define the evaluation of spin networks and construct the 
whole picture behind the paper. It is in the same spirit of this section that all 
of the future sections follow. 

2.1 500(2, 1) and its principal unitary representations 

Recall 3-dimensional Minkowski space-time to be M 3 with its Lorentzian metric 
given by 

[x, y] = x y - xiyi - x 2 y 2 

where x = (xo,X\,X2),y = (yo,yi,V2) £ and xq and yo are the time compo- 
nents. 

Consider the following subsets of Minkowski space: 

1) The Cone C: This is the set of points of Minkowski space which satisfy 
[x, x] = 0, that is, the set of vectors which are at a distance zero from the origin. 
These are called the lightlike vectors as this set describes particles which travel 
at the speed of light. 
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The upper cone C + : This is the set of points in C for which the time 
component x > 0. 

2) The two-dimensional hyperbolic space H 2 : This is the set of points for 
which [x, x] — 1 and xq > 0. 

Consider SO (2, 1), the group of unimodular linear transformations of Minkowski 
space-time which preserves the Lorentzian metric. The subgroup of 5*0(2, 1) 
which preserves C + is denoted 5*00(2, 1). This group is connected and locally 
compact. 

In fact the action of SOq(2, 1) in M 3 gives an action of 500(2, 1) in the 
hyperbolic space H 2 , and the subgroup of 50o(2, 1) which leaves the vector 
(1,0,0) fixed, is isomorphic to 50(2). It is a maximal compact subgroup and 
therefore H 2 can equivalently be seen as a homogeneous space given by the 
quotient 5O (2, l)/50(2). 

The representations of 50o(2, 1) are realised in the following way. First of 
all, we introduce spherical coordinates on the hyperbolic space H 2 . That is, for 
x = (xo,Xi,x 2 ) G H 2 its coordinates are given by 



x = coshr 
x\ = sinh r sin 6 

X2 = sinh r cos 6 (5) 



In these coordinates the Lorentzian wave operator given by 



d 2 d 2 d 2 

+ 77^ + 773 ( 6 ) 



Oxq dx\ dx\ 
when restricted to the hyperbolic space H 2 is given by 

d d d 2 

A + = (sinhr) -1 — sinhr— — h (sinhr)" 2 — — - (7) 
Or or o9 2 

Consider the space of functions inside C + . The group 50o(2, 1) acts on this 
space and the operator of equation (7) commutes with this action. Then if there 
is a function / for which Af(x) = then Af(g~ 1 x) = 0, for g G 5 (2, 1). A 
function which satisfies A/(x) = is called harmonic. Consider the space of har- 
monic functions defined inside C + of degree a G C. The equality T a (g)f(x) = 
f(g~ l x) defines a representation of 50o(2, 1) in the space of harmonic func- 
tions. There are the ones called principal representations which are the ones 
we are interested in in this work. These representations are labelled as (ip — |), 
where p is a real parameter. Its importance for the physics of this work is that 
tetrahedra embedded in 3-dimensional Minkowski space for which edges are la- 
belled by principal unitary representations are thought of as corresponding to 
spacelike tetrahedra, that is all their edge vectors are spacelike. 
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2.2 The kernel and its asymptotics 



As in jHj we approach to the calculation of the kernel which is the main ingredient 
for the evaluation of the spin networks, and most important to our case of the 
evaluation of the tetrahedron. 

Given a function h G L 2 (H 2 ) we have the p irreducible component 

f p (0= [ h(x)[x,^dx (8) 

where £ is a lightlike vector, that is £ G C + . According to j!4j . the inverse 
transform is given by an integral of the type 

h(x) = / dn p / f p {£)\z,€\-t-+>d£ (9) 
Jo Js 1 

where S 1 is a contour on C + intersecting every generatrix of the cone at one 
point. More specifically for our purpose it is taken to be the unit circle S 1 . 
Composing these two transforms without integrating over p gives the projection 
operator onto the p irreducible component h p given by 

h P (x) = ^-[ ! [x,^- ip [y,^ +ip h(y)d^dy (10) 
4«vr J H 2 J s i 

where we consider the kernel function 

K p (x,y) = [ [x,^- ip [y,^ +ip dC (11) 
Js 1 

where x, y G H 2 , and then 

h P (x) = ^~f H2 k p& y)Hy)dy (12) 

We now consider the kernel and express it in a spherical coordinate system. For 
this to be done, we shift the points x, y G H 2 (which always can be done) to 
the following points 

x = (cosh r, 0, sinhr) y — (1,0,0) 

by abuse of notation we still call them x and y. Moreover coshr = [x,y] is 
the hyperbolic distance between x, y G H 2 . The point £ G 5* 1 has coordinates 
(1, sin#, cos6*). Hence we write 

1 f 2?T 

K p (x,y) = — / (coshr — cos(a) sinh r)~ 2 ~ ip da (13) 
2tt Jo 

It turns out that our kernel is indeed the well known zonal spherical functions 
for representations of the group SOq(2, 1), known also as Legendre functions 
*B T (cosh r) where r = — 1/2 — ip. Some of its properties are described as follows. 
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This kernel is an eigenfunction of the Laplace operator on H 2 . That is, 
A + K p (x,y) = —(1/4 + p 2 )K p (x,y). Because of this property, the function 
K p (x,y) is called a spherical harmonic. 

Now, given a spin network with edges labelled by real numbers M + cor- 
responding to principal unitary irreducible representations, we evaluate it as 
follows. To each vertex we associate a variable x G H 2 . To each edge labelled 
by p, we associate the kernel K p (x,y) where the variables x and y are the ones 
associated to its correspondent vertices. 

The evaluation is then given by the integral 



where, as in we have regularized the evaluation by removing the integration 
over one of the variables in H + at one of the vertices. This variable is then 
fixed. 

From this definition we can attempt to evaluate some spin networks. It 
follows easily for example that the evaluation for our face amplitude 



Some other examples such as our edge and vertex amplitudes result in integrals 
which are difficult to evaluate directly. So we overcome this problem by consid- 
ering the asymptotics of one of these evaluations, which is of our main concern. 
We then study the asymptotics of the tetrahedron network. 

Given our tetrahedron network labelled by principal unitary irreducible rep- 
resentations Pi2,Pi3,Pi4,P23,P24,P34 as in figure 1. 




(14) 




gives 



K p (x, x) — 1 



X 




Figure 1: Tetrahedron network 



Naively we evaluate it as the following integral 
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T 6 = / dx 2 dx 3 dx4K Pl2 (x 1 ,X2)K Pri (x 1 ,x 3 )K Pl4 (x 1 ,x 4 ) 

JH 2 xH 2 xH 2 

x K p23 (:c 2 , a*)-^ (z 2 , x 4 )ir p34 (x 3 , x 4 ) (15) 

Note that this is an integral involving other integrals, that is, the kernels which 
are also given by integrals. We therefore have to know whether our tetrahedron 
evaluation converges. We do not have a way to prove that our tetrahedron 
integral is finite but we give some evidence about it. First of all the problem 
is very similar to the 4-dimensional Barrett-Crane model of evaluating the lOj 
symbol, and this is being proved to be finite. 

We can prove the convergence of the theta symbol. This latter is given as 
follows: We have that our Lorentzian kernel when the distance r between the 
points x and y goes to infinity and p ^ 0, [13J, is asymptotic to 



K p (x,y)~A 7W — (16) 



1 

(2cosli"r) 1 /2 

where A is a factor which oscillates. 

Then for very large r we have that our kernel is asymptotic to ~^ e ~ r ■ With 
this asymptotics our theta symbol evaluation is then given by 




N ww e ' ir (17) 

which is finite. 

We then have a finite value of our theta symbol and this makes the conver- 
gence of our tetrahedron more promising. More evidence is given is section 4 
when dealing with our Newtonian quantum gravity model; the Lorentzian and 
Newtonian models are related by a limit as explained in section 4. 

Supposing then that our tetrahedron evaluation converges, we study the 
asymptotics of this integral. 

As a warm-up example let us compute the asymptotics of our kernel K p (x, y) 
first. 

Recall that the general method of stationary phase states that an n-dimensional 
integral of the form 

I{x) = J g( x )e iaf{x) dx (18) 
has an asymptotic expansion when a — > oo given by 

-2ttW 2 ^ , , inff ^ e ^sgn{H)/A 
x\df(x)=0 



/9-7r\«/2 p insgn{H)/A 

'<*>-(£) E sM^ W ^m F + °(«-" /2 - 1 > (») 
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where H(x) is the Hessian matrix for f(x), and sgn(H) is its signature, i.e the 
number of positive eigenvalues minus the number of negative eigenvalues. 

Now observe that for the case of our kernel we have g(x) = (coshr — 
cos (a;) sinhr)^ and e ik ^ = e v(-Mcoshr-coe( a! ) B iiihr)) ? g0 that = _ hi(coshr- 
cos(x) sinhr). Then 

/ sinhrsin(x) 

coshr — cos(x) sinhr 

and so the stationary points of are given by x = and x = ir if r 7^ 0. 

Then following the stationary phase formula, the asymptotics of K p (x,y) as 
p — > 00 is given up to a 1/27T factor by 

p ' Vpsinhr/ L (cosh r — sinhr) v (coshr + sinh r) tp J 

Using trigonometric identities it is easy to see that our above formula reduces 
to 

K P (x, y) ~ ( — ^— ) 5 U^/ 4 ) + e-^^-^l (22) 
Vpsinhr/ L J 

which can also be seen as 

Kp{x ' V) ~ (^7) 3 C ° s(pr _ ^ (23) 
2.3 The tetrahedron network 

On the assumption that our tetrahedron network is finite, we now proceed to 
compute the asymptotics of formula (2). For that purpose we replace the labels 
Pij by apij to study the asymptotics when a — > 00. We then have 

T 6 = / TT/ (coshry -cos(%)sinhr ii )-5+ i ^^. j rf Xj (24) 

which we now write as 

x J e iaE i<3 ^ln(cosh^-cos(^) S inh^) d Q. jdx . ( 2 5) 



The exponent 5 = Yli<j Pij ln(coshry — cos(^) sinhr^) in the equation above 
is called the action from which we have to compute the stationary points in 
order to obtain the asymptotics of our tetrahedron network. 
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Solutions for which our action is stationary should be classified as either 
degenerate ones or non-degenerate. Degenerate solutions correspond to a de- 
generate tetrahedron, that is, one for which any or all of its vertices coincide, 
Xi = Xj for any % 7^ j, or equivalently when some or all r^- = 0. 

Non-degenerate solutions correspond to tetrahedra for which 7^ for all 

We then start our study of all these possible solutions and follow a procedure, 
similar to the one used in [I] 



2.3.1 Non-degenerate solutions 

We have in this case that all of the points X{ G H 2 are different, that is 7^ 
for all 

We start varying the action 

S = J2i<j Pi j In (cosh rj j — cos(9ij) sinhr^) with respect to and Xj. We 
use a Lagrange multiplier Xj for the constraint we have [a^a^-] = 1. 

A first set of equations are given by varying our action with respect to 
that is 



OS 

from which it follows 



<9% 



Pij sin(^y) sinhr 



'j 



coshrjj — cos(%) sinhrjj 



(26) 



(27) 



Then the above equation, for the case for which r$,- 7^ for all i 7^ j, is satisfied 



when 9ij = or = n 

To keep track of the possible choices of 9ij = or 9ij = tc we introduce a 
variable ey = — 1 if 0^- = or = 1 if 0^- = 71. 

Our action can then be written S = J2i<j e ijPij r ij so that a secon d set of 
equations are given by varying our action with respect to 22, £3, ^4 that are the 
variables of our integrand. We take care of our constraint mentioned before, 
that is 

OS 

— = XjXj (28) 

where j = 2, .., 4. 

From this it is easy to see that we arrive at the following equations 



e 12Pl22 ; l . C23P23 a; 3 . e24P24^4 ^ 

~^~~r 1- ~^—r 1- ~^—r = ^2X2 

smn ri2 sinn r23 smh 

ei3Pl3^1 , e23P23^2 , C34P34a;4 . 

h h ~^~T = A 3^3 

smh r 13 smh r 2 3 smh r 34 
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ei 4 pi4^i 



+ 



e 2 4P24^2 



+ 



e 3 4P34X 3 

sinh r 34 



A4X4 



(29) 



sinh ri4 



sinh r 2 4 



where we recall that all of our vectors Xi are timelike and future directed, and 
our e are the result of the possible election of the 9 stationary points at or n. 

We now take the wedge product with xi to the i equation(? = 2, ..4), 2 where 
we take a unit vector in the direction x^ A Xj. So we have the following set 
of equations, 



If we sum the three equations above, we get a fourth equation given by 



If we call Vij = eijPijVij, then the equations (28) and (29), imply that the are 
the edge vectors of a tetrahedron with edge lengths pjj. Each equation describes 
a face of a tetrahedron. We have that the normals to each of these faces are 
timelike vectors, therefore our tetrahedron is spacelike which implies then that 
the usual triangle inequalities should be satisfied for the lengths of its edges. 
Each of our timelike normal vectors could be pointing inwards or outwards, 
but always future directed. Let rii be the outward normal pointing vectors to 
the faces of the tetrahedron. Then we have that where ai = ±1. 

Recall that our point x x is fixed at (1, 0, 0). 

All the possible configurations of a spacelike tetrahedron are divided in two 
possible configurations only: 

1) One normal vector, say ni, is future directed, while the remaining three 
normal vectors Xi for i — 2, 3, 4 are past directed. 

2) Two normal vectors are future directed and two normal vectors are past 
directed. 

We could also have one normal vector, say n 1; is past directed, while the 
remaining three normal vectors Xi for % = 2,3,4 are future directed. But this 
latter configuration is equivalent to the first one as we need only to change all 
the signs. 

We then have that the possible values of the e's are: 

For configuration 1), we have that x± = ni,x 2 = — n 2 ,x 3 = —713,2:4 = — n 4 , 
therefore e 12 = -1, 613 = -1, £i4 = -1, £23 = 1, <^4 = 1, e 34 = 1. 

For configuration 2), we have that X\ = rii,x 2 = n 2 ,xs = — n^^x^ = — n 4 , 
therefore e 12 = 1, ei 3 = -1, e u = -1, e 2 3 = -1, e 24 = -1, e 34 = 1. 

2 Recall that, (xj A Xi) a = T] ae e e b c x b jX'i, where e a b c is the totally antisymmetric tensor with 
£123 = 1 and r/ a b = rf h is the — V + Lorentzian metric. 



- ei 2 pi 2 -ui 2 + e 23 p 23 i; 2 3 + e 24 p 24 i; 2 4 = 

-ei3Pl3^13 - C23P23^23 + e 34 P34^34 = 
-ei 4 pi 4 f 14 - e 24 p 2 4^24 - e 34 p34^34 = 



(30) 



ei2Pl2fl2 + ei3Pi 3 Wi 3 + ei 4 pi 4 Vi4 = 



(31) 
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Other equivalent configurations are given by swapping all the signs of the 
above configurations. 

The Hessian is a 11 x fl matrix with components 9 ®.q x . , dx . d S g .. , 90 89 ■■ ■ 
When evaluating the Hessian at the stationary points, many of the entries will 
have zeros and it is seen that its determinant equals the determinant of the 
5x5 matrix with components 9 °.q x . evaluated at the stationary points, times 
some factors. The remaining factors are given by (pi2---P34)/(Yli<j coshr^- — 
cos(0y) sinhr^) - ^) which denominator cancels in the evaluation of the asymp- 
totic expansion. 

Our integral is of dimension 11 which then implies that there is an extra 
factor of (2n/a) u / 2 

The contribution we will get will be an oscillating function as 

T 6 ~A cos(J2 Pa r ij + k ( 32 ) 
k<i 

There is a factor A containing information about the determinant of our 
Hessian at each stationary point. 

The information contained in the factor A is related to the volume of the 
tetrahedron. This specific calculation can be done numerically for some exam- 
ples, but a rigorous proof must be done in a different way. 



2.3.2 Degenerate solutions 

We now proceed to the study of the contribution which comes from degenerate 
solutions. Let us deal with the case in which all our vectors x's are parallel. We 
have noted before that our kernel was not very promising when dealing with 
this case. However, we change the way our kernel looks so that we may deal 
with our degenerate case. For this purpose we introduce new coordinates as 
follows: given two points Xi and Xj in H 2 there is a geodesic j(t) which joins 
the points Xi and Xj. We define the vector ^ to be the initial velocity of our 
geodesic. Then we have | ^ \— rij is the distance between Xi and Xj. Recall 
our kernel 

/•27T 

K p (xi,Xj)= / (coshr — cos(a) sinhr)~ * +ip da (33) 
Jo 

Let yij be unit tangent vectors, that is, G S 1 , so that cos(a) = (£^/ | ^ \)-yij. 
As the function sinh r/r converges to 1 as r — > we change sinhr for r, which 
is just given by | ^ \. Then our kernel can be written as 

K p (xi,Xj)= {[xi,Xj\-£ij-yij)~^ +%p dy i: j (34) 
Our action S is given by 
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Now if we first vary the x's and then and then make the x's to be parallel, that 
is &j — > 0, we obtain the equations 

- P12Z/12 + P23Z/23 + P242/24 = 
-P132/13 - P23Z/23 + P342/34 = 

-PUVU - P242/24 - P341/34 = (36) 

If we sum these three equation we get a fourth one given by 

P122/12 + P132/13 + PuVu = (37) 

Each of our equations (33) and (34) describe triangle faces, and the vectors yij 
are two dimensional, as they are vectors in S 1 . Then the equations describe 
the faces of a degenerate tetrahedron. This means that our tetrahedron lies 
completely in a two-dimensional plane. Once the lengths of edges, are specified 
the geometry of the tetrahedron is completely specified. Therefore the above 
equations are not satisfied generically, which implies that we cannot expect a 
contribution from degenerate solutions of our tetrahedron evaluation. This is 
analogous to the Ponzano-Regge result, in which they conclude that the only 
contribution to the 6j-symbol comes from non-degenerate solutions. We have 
a similar argument for our Euclidean case of section 4, so that once more we 
agree with Ponzano-Regge. 

Conclusion. As a conlusion to this entire section we conclude that the only 
contribution we get is from non-degenerate tetrahedra which in accordance to 
our asymptotic formula (32), and to the result of the asymptotics of the squared 
6j-symbol of jTHj we have the asymptotic formula 

T 6 ~ -7= COs(^2 PijTij + k^j (38) 
V V k<l 

where V is the volume of our tetrahedron. 

3 Euclidean spin networks 

With the same technique that we developed for the Lorentzian case, we can 
deal with the Euclidean spin networks. Equation (3) is the partition function 
for (2+l)-dimensional Euclidean quantum gravity as a spin foam model, and it 
is equivalent to the Ponzano-Regge model. 
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3.1 SO(3) and its unitary representations 

We now consider the 3-dimensional Euclidean space-time, which is M 3 with the 
usual positive definite metric given by 

[x, y] = x y + xij/i + x 2 y 2 

Consider the two dimensional sphere 5 2 , which is the subset of M 3 given by the 
set of points for which [x,x] = 1. 

SO (3) is the group of unimodular linear transformations which preserve the 
positive definite metric. It is compact and connected. The stabilizer of the unit 
vector (1,0,0) is isomorphic to 5*0(2). Therefore our unit sphere S 2 can be 
seen as the homogeneous space 50(3)/50(2). 

The representations of 5*0(3) are realised as follows. Introduce spherical co- 
ordinates on the unit sphere S 2 . That is, for x = (x , x 1: x 2 ) £ S 2 its coordinates 
are given by 



x = cos 9 
x\ = sin 9 sin ip 

x 2 = sin 9 cos (p (39) 



In these coordinates the Euclidean Laplace operator 



d 2 d 2 d 2 

+ 7^ + 7^2 (4°) 



dxl ® x l ® X 2 
when restricted to the unit sphere S 2 is given by 

A = (sinh^)- 1 — sin^— + (smlir)" 2 — (41) 

Consider the space of functions on M 3 . The group 50(3) acts on this space 
and the operator of equation (40) commutes with this action. Then if there 
is a function / for which Af(x) = then Af(g~ 1 x) = 0, for g G 5(3). A 
function which satisfies A/(x) = is called harmonic. Consider the space of 
homogeneous harmonic polynomials of degree £ in the variables xo,xi,x 2 . The 
equality T e (g)f(x) = f{g~ 1 x) defines a representation of 50(3) in the space of 
homogeneous harmonic polynomials. These are unitary representations and are 
labelled by positive integers £ = 0,1, 2, .... 



3.2 The kernel 

The question now is, how does the kernel look like? We recall that our Lorentzian 
kernel is given by a zonal spherical function of the corresponding group 50o(2, 1), 
more specifically a Legendre function *B_i/ 2+ j P (coshr). 
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It is then obvious to think that for the Euclidean case the kernel is given by 
a zonal spherical function of 5*0(3). This kernel is then given by the integral 



1 f 27T 

K £{x,y) = 7T / ( cosd 



+ i cos <f sin 9) dip 



(42) 



where x,y G 5 2 . 

The above kernel is the well known Legendre polynomial Pe(cos8) where 
now £ denotes an irreducible representation of SO (3). 

It happens that the kernel (42) has the following asymptotics 



3.3 The tetrahedron network 

We define the evaluation of Euclidean spin networks in the same way we did 
it for the Lorentzian case but now using our kernel given by the the Legendre 
polynomial Pi(cos9). We can also study the asymptotics of the tetrahedron 
network for this case. The way to proceed is as we did it for the Lorentzian and 
Newtonian case, the only difference being that we now have discrete irreducible 
representations L 

Given our tetrahedron network labelled by unitary irreducible representa- 
tions of 5*0(3), as in figure 2. 



we have its evaluation given by the integral (4). 

We can as well divide the study of our solutions in non-degenerate and 
degenerate ones. According to Ponzano-Regge ^2], there is no contribution 
from degenerate configurations for this Euclidean case. We should be able to 
visualise this with the method we follow. 

3.3.1 Non-degenerate solutions 

The tetrahedron network now can be defined by the following integral 



/ 2 \ 1 / 2 / 1 7T\ 



(43) 




Figure 2: Euclidean tetrahedron network 
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T 6 = I T\f gaEi^^Mcos^+icos^^sinh^)^.^ ^ 

-As 2 ) 4 7<j J l°M 6 

Notice that in this case we have an integration carried out over 4 copies of 
the sphere S 2 , one for each point of our tetrahedron. This is done because for 
this case we have a compact space S 2 , and then there is no divergence when 
integrating over these 4 copies. There is also a very slightly difference here 
that is interesting to notice. Our Euclidean tetrahedron integral has a complex 
integrand. From first instance it appears that the asymptotics of our integral 
has to be treated by the steepest descent method. However if we take that 
our action is given by S = X)i<j ^ij ln(cos0jj + % cos(tpij) sinh#y). we first of all 
notice that if we vary the action with respect to the angles </?'s we find that 

dS _ -iiy sin ipy sin 9 i:j 



difij cos 9ij + i cos ifij sin 

Its stationary points are given by ifiij = and ipij = n. Introducing a variable 
€ij = 1 if (fiij = 0, or €ij = —1 if ip^ = it, our action can now be written as 
S = i Yli<j e ij^ij@ij- Then as an i term has appeared before our action, it means 
that now we can deal with the integral (44) by using the usual stationary phase 
method. 3 

Then, we can vary our action with respect to the four rr's variables, paying 
attention to the Lagrange multiplier obtaining the following equations 



€12^12^2 . €13^13X3 61^14X4 

■ 2 + —3 + • a = X ^ 

sin & u sm &i 3 sin 14 

£12^12^2 , £23^23^3 , £24^24^4 x 

— a — + — a — + —7 — = A 2 x 2 
sm 6*12 sm #23 sm # 24 

£13^13^1 , £23^23^2 , e 34 £ 34 X 4 x 
• 2 + —3 + • a = A 3^3 

sm t7 13 sm t^23 sin #34 

euiuXi 624^24^2 . ^34^34^3 x / Ar ,\ 

■ + • Q + • Q = ^4^4 46 
sm &i4 sm t>24 sm 6/34 

We can then take the Euclidean wedge product with X; L to the i equation 
(i=l,..4), where we take a unit vector Vij in the direction X{ A Xj. Then we 
get the following equations 



£12^12^12 + £13^13^13 + £14^14^14 = 
-£12^12^12 + £23^23^23 + £24^24^24 = 
"£13^13^13 - £23^23^23 + e3 4 £ 34 W 34 = 

-eidi A V XA - e 24 ^24^24 - £34^34^34 = (47) 



3 We will notice that this cannot be done when dealing with 'degenerate configurations'. 
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If we call Vij = 6ijPijVij, then the equations (47) imply that the V^- are the edge 
vectors of a tetrahedron with edge lengths £ij. Each equation describes a face 
of a tetrahedron. Each of our Xi normal vectors could be pointing inwards or 
outwards. Let n« be the outward normal pointing vectors to the faces of the 
tetrahedron. Then we have that where cij = ±1. Recall that our 

point X\ is fixed at (1, 0, 0). 

All the possible solution configurations are being classified in [3]. Moreover, 
all of our representations are labelled by integers, which means that the sum of 
the labels at each vertex is an integer, meaning that all stationary points either 
give the same contribution or give the complex conjugate one. Then given a 
solution the corresponding asymptotics is given by 

T Q ~A cosQT(^ + ^)% + £;£) (48) 

k<l 

where again A contains information about the determinant of the Hessian, and 
k is the signature of it at a stationary point. 

The information contained in the factor A must again be related to the 
volume of the tetrahedron. This specific calculation can be done numerically 
for some examples, but a rigorous proof must be done in a different way. 

Degenerate configurations. Degenerate configurations are treated in the 
same terms we dealt with the Lorentzian degenerate configurations. All the 
calculations are similar, and the conclusion to which we arrive is that no con- 
tribution is really expected from these configurations. 

4 Newtonian spin networks 

Following a similar fashion to our previous Lorentzian and Euclidean case, we 
deal now with which we call Newtonian. This name is suggested from 

the certain limits which will be understood in the present section. 

First of all we can just make notice that our Lorentzian and Euclidean kernels 
have a similar asymptotics to the 0th order Bessel function Jo(z). The 0th order 
Bessel function has an integral representation given by 

J Q ( Z ) = — [ e lzcos{e) d6 (49) 
2tt Jo 

which asymptotics is given by 




This 0th order Bessel function is a zonal spherical function for representations 
of the group I SO (2) which is the group of motions of the plane, that is, the 
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group of nonhomogeneous linear transformations of the plane which preserves 
distances between points and the orientation. 

Mathematically ISO{2) can be obtained by a limit procedure from SOq(2, 1) 
|14j . From this limit the principal unitary irreducible representations of ISO(2) 
are again given by real numbers R + and then our previous kernel K p is given by 
the Oth order Bessel function Jo(pr) where now our r distance between points 
is given by the usual flat 2-dimensional metric on a plane. 

Physically this limit corresponds to taking the limit in which the speed of 
light goes to infinity (c — > oo) . Then it is natural to think of our Bessel function 
as a kernel, which gives us a way to compute spin networks for a Newtonian 
quantum gravity model. 



4.1 ISO (2) and its irreducible representations 

Physically we can understand the Newtonian model by consider again Minkowski 
space-time R 3 with Lorentzian metric, and then taking the limit in which the 
speed of light goes to infinity. Then our metric turns into a degenerate metric 
with signature (0,+,+). From this limit our hyperbolic space H 2 becomes a 
flat plane isomorphic to R 2 , and our Lorentz group SO (2, 1) become the group 
ISO(2). 

ISO(2) is the group of motions of the plane R 2 . It is the group of nonho- 
mogeneous linear transformations of the plane, and preserves distances between 
points, as well as the orientation of the plane. Every element of ISO{2) is a 
composition of a rotation about a point generated by an element of the group 
50(2), and a translation along a vector. We have that our flat plane M 2 is seen 
as the homogeneous space ISO{2) / SO{2). 

The Laplace operator in R 2 is given by 



dx\ dx\ 

which in spherical coordinates is given by 



A = 7^ + 7w2 (51) 



a 92 19 1 92 

^=Q^ + -Q- r +-2Qf2 (52) 

The irreducible representations of IS 0(2) can be realised by considering the 
space of functions on R 2 satisfying the equation Af(x) = —R 2 f(x). The equa- 
tion 

(T»(g)f)(x) = f(g- 1 x) (53) 

for p = iR, defines a representation of the group I SO (2) in the space of functions 
on R 2 , which is irreducible for R ^ 0, and unitary when R 6 R. 

Moreover, the Bessel functions are eigenf unctions of the Laplace operator 
on the plane R 2 . 
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4.2 Spin network evaluations 

Before we study the tetrahedron network we give some examples of the evalua- 
tion of spin networks in the Newtonian quantum gravity model. 

A Newtonian quantum gravity model in 2+1-dimensions evaluates ampli- 
tudes of spin networks in the 2-dimensional plane M 2 , and the evaluation of 
these are given by using a propagator in this plane which is given by the 0-th 
order Bessel function. We follow the usual recipe for evaluating the amplitudes 
of spin network graphs which for our case is given by 



A(graph) = I / dx v \\ J (p e r(x,y)) (54) 

p e denotes a continuous representation of the group I SO (2) which labels an 
edge, and r(x,y) is the usual distance between points x and y in a two dimen- 
sional plane. 

Examples of the evaluations are given by 




pi 

5{p 2 -pi) 




pi 

1 

A 

where A is the area of the triangle formed by edges of length side pi, P2, P3- If 
not triangle is formed, the evaluation of the theta graph is zero. This shows 
once more that a convergence of the Lorentzian tetrahedron network is likely 
to be expected, since Lorentzian quantum gravity is related to the Newtonian 
case, and we can see that for the Newtonian quantum gravity model the exact 
evaluation of the theta graph exists. 

4.3 The tetrahedron network 

We can now define a way to compute a tetrahedron network for this Newtonian 
case, and study its asymptotics as we did with the Lorentzian case. 

Moreover we can attempt to define a spin foam model of Newtonian quantum 
gravity accordingly to the weight amplitudes given to our vertex, edge and face 
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spin networks. Particularly we have that a Newtonian tetrahedron network has 
an evaluation given by the integral 



T 6 = / dx 2 dx 3 dx 4 J (p 1 2r(x 1 ,x 2 ))Jo(pi3r(x 1 ,x 3 ))J (p u r(x 1 ,x i )) 

</(R 2 )3 

X J (P23?" (X2 ,X 3 ))J (p 2 AT (x 2 , X 4 ) ) Jq (p34?~ (x 3 , X 4 )) (55) 

Pij denotes again a continuous representation of the group I SO (2), and r(xi, Xj) 
is the euclidean distance between points Xi and Xj in a two dimensional plane. 
We simplify our notation by denoting r\j for the distance between our points 

Xi 5 Xj . 

We can then substitute our integral representations of our Bessel functions 
and get 

T e = [ TT I e i ^ pi ^ cos ^ ) d9 ij dx j (56) 

i(R2)3 f<) J[0,2nf 

For this case our points Xj, for j = 1,,4 now live on an embedded flat 2- 
dimensional plane at our 3-dimensional space-time at time 1, where our metric 
turned into a degenerate one but when restricted to our 2-dimensional plane, it 
is just the usual flat Riemannian metric on that plane. 

We proceed similarly as in the Lorentzian and Euclidean cases. 



4.3.1 Non-degenerate solutions 

The exponent S = X^<j Pn r n cos (%) m the equation above is called the action 
from which we have to compute the stationary points in order to obtain the 
asymptotics of our Newtonian tetrahedron network. 

It is again a similar procedure as the one we followed for the Lorentzian 
and Euclidean case, where we find the asymptotic points by varying our action 
with respect to the variables Xj, and For the 9 variables we just make the 
variation to be zero. As for the x variables, we can think of them as vectors 
with coordinates x = (1,X), where X is a vector in the M 2 plane. The action 
can be considered as a function of the 9 variables and of these two dimensional 
vectors X. 

We now look at the stationary points. 

The equations we obtain are analogous to the ones we had for the Lorentzian 
case and are given by a much simpler formula 

dS 

— = - sin 9 tj (57) 

which stationary points are given by 9ij = and % = tt. 

Then our action can be written as S = X^<j e ijPij r ij where is given just 
in terms of the distance in the plane 1R 2 by = ((Xi — Xj) ■ (Xi — Xj)) 1 / 2 . 
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Varying our action with respect to the variables X 2 , X 3 , X 4 , gives the equa- 
tions 



(X 1 - X 2 ) 

. e 23P23 (X 2 - X 3 ) 

. ^2Ap2i{X 2 — X4J 

r\2 r 23 r 24 

e 13Pl3 {X 1 -X 3 ) £23P23(X 2 — X 3 ) e 34 p 34 (X 3 — X4) 

ri3 r 23 r 34 

(*i - x,) 

. ^2AP2a{X 2 — X4) 6 34 p 34 (X 3 — X 4 ) 
ru r 24 ?"34 

We note that a solution to the equations (58) is given by the appropiate signs 
given to the e's. Given such solution, interchanging all the signs of the e's gives 
the complex contribution. Then, we can expect an oscilatory contribution in 
the cases for which there is a solution to the equations (58). Moreover, as in this 
case we have a degenerate metric, a degenerate tetrahedron is already implicit 
in the theory. A better understanding of all the relations between the Netonian 
theory and Lorentzian or Euclidean quantum theory is an interesting problem 
to work in. This also applies to other dimensions as mentioned already in pQ 
for the 4-dimensional case. 



= 
= 

= (58) 



5 Conclusions 

We have observed that all of our kernel functions of the Lorentzian, Euclidean 
and Newtonian models are given by zonal spherical function of a representation 
of the respective group. In particular we think of the Newtonian the limit 

case of Lorentzian quantum gravity in which the speed of light tends to infinity. 
This gives the idea that we have a picture of a unified theory of quantum gravity 
and spin networks. 

We could generalise all the ideas of the paper in order to deal with spin 
networks in any dimension. We therefore have that Lorentzian spin networks 
for n-dimensional quantum gravity whose group is SO (n — 1,1) should be 
evaluated by using a kernel given by a zonal spherical function of this group. 
That is, 

Kj;(r) = V f' / (cosh r - cos ^sinhr) CT sin"" 3 QdO (59) 

where a = —p + ip. p is related to the dimension of the space-time as p = 
(n-2)/2. 

Euclidean spin networks for n-dimensional quantum gravity whose group is 
SO(n), should then be evaluated by using a zonal spherical function of this 
group. That is, 
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jy re— 1 > 

iff (0) = 1 2 n _ 7 / (cos 9 — % cos <£> sin sin n - 3 pefy> (60) 
y/nT(— ) Jo 

And finally Newtonian spin networks for n- dimensional quantum gravity whose 
group is ISO(n), should then be evaluated by a zonal spherical function of this 
group. That is, 

TV re— 1 \ p% 

K?(r) = J_ 2 J 2 / e^ ose sm n ~ 3 ede (61) 
v 71 " 1 {—) Jo 

The above functions are Legendre *B"(coshr) functions, Legendre polynomi- 
als P™(cos6 l ), and Bessel functions J n (pr) respectively, of the corresponding 
dimensions. 

We can see that for 3-dimensions we have again the kernels we studied in 
this paper. 

For the 4-dimensional case, we have a = —1 + ip, that is, a principal unitary 
series representation of the Lorentz group 5*00(3, 1), then our Lorentzian kernel 
reduces to 

«?«-3£ < 62 » 

which is the Barrett-Crane one. 

Similarly for the Euclidean and Newtonian kernels, we get 

Kf«f) = p^± 1V (63) 

K»At) = ^ (64) 

The way all these kernels are related is given in [Hj. In particular for our 
3-dimensional case we have that 

lim K L {rjt) = J (r) (65) 
lim Kf{B/l) = J (r) (66) 

£— >oo 

This is interpreted as saying that when dealing with the degenerate configu- 
rations in the Lorentzian and Euclidean spin network evaluations, these can be 
related to the Newtonian case. A better understanding of what all this implies 
to the models is an interesting problem. In particular what we can now about 
the physics of these models is important. We want to stress once more the 
interest in the Newtonian model. In the end, it is related to the Lorentzian and 
Euclidean, and it is simpler. It is then interesting to study Newtonian quantum 
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gravity as itself and try to find what it could be telling us about the physics of 
quantum gravity. 
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